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We report on numerical lattice QCD calculations of some of the low moments of the nucleon 
structure functions. The calculations are carried out with gauge configurations generated by the 
RBC and UKQCD collaborations with (2+l)-flavors of dynamical domain wall fermions and the 
Iwasaki gauge action (/3 — 2.13). The inverse lattice spacing is a~^ = 1.73 GeV, and two spatial 
volumes of (2.7 fm)^ and (1.8 fm)^ are used. The up and down quark masses are varied so the pion 
mass lies between 0.33 and 0.67 GeV while the strange mass is about 12 % heavier than the physical 
one. The structure function moments we present include fully non-perturbatively renormalized iso- 
vector quark momentum fraction, {x)u-d, helicity fraction, (a:)Ati-Ad, and transversity, {l)su-Sd, as 
well as an unrenormalized twist-3 coefficient, d\. The ratio of the momentum to helicity fractions, 
{x)u-d/ (x) Au-Ad, does not show dependence on the light quark mass and agrees well with the value 
obtained from experiment. Their respective absolute values, fully renormalized, show interesting 
trends toward their respective experimental values at the lightest quark mass. A prediction for the 
transversity, 0.7 < {l}su-Sd < 1.1, in the MS scheme at 2 GeV is obtained. The twist-3 coefficient, 
di, though yet to be renormalized, supports the perturbative Wandzura-Wilczek relation. 

PACS numbers: ll.15.Ha, ll.30.Rd, 12.38.Aw, 12.38.-t 12.38.Gc 



I. INTRODUCTION 

We report numerical lattice quantum chromodynam- 
ics (QCD) calculations of some low moments of nucleon 
structure functions using the lattice gauge ensembles [1] 
jointly generated by the RIKEN-BNL-Columbia (RBC) 
and UKQCD Collaborations with "2-hl" flavors of dy- 
namical domain- wall fermions (DWF) PHI]- Recently, 
there has been an increased interest in lattice calcula- 
tions of these moments (see [SHZ] for recent reviews) . 

The structure functions are measured in deep inelas- 
tic scattering of electrons off a nucleon , the cross 
section of which is factorized in terms of leptonic and 



hadronic tensors, cx la/^W"^ . Since the electron leptonic 
tensor, l^is, is known, the cross section provides us with 
structure information about the target nucleon through 
the hadronic tensor, 

= i J d*xe"^''{N\Ty{x)j'^{0)]\N). (1) 

Here q denotes the spacelike four-momentum transferred 
to the nucleon from the electron through a virtual pho- 
ton. The hadronic tensor is decomposed into symmetric 
unpolarized and antisymmetric polarized parts, W"^ — 
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pa 



(2) 



with kinematic variables defined as Pa the nucleon mo- 
mentum, Sa the spin normahzed with the nucleon mass 



q- P, X = Q^/2v, and 



The unpolarized structure functions are Fi{x,Q ) and 



92{x,Q )] 5 92(x,Q , 



(3) 
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F2{x,Q'^), and the polarized, gi{x,Q'^) and g2{x,Q^). 
Their moments are described in terms of Wilson's op- 
erator product expansion: 
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where the Wilson coefficients, ci, C2, ei, and 62, are per- 
turbatively known. The moments, {x^)q{^), (x")ai3(m) 
and (/i) are calculable on the lattice as forward nucleon 
matrix elements of certain local and gauge-invariant op- 
erators. 

In addition, the tensor charge. 



M 

1)5<?(m) = o/c p ?r^-T(^,'5|gicra/375g|-P,S'), 



(8) 



which probes the transverse spin structure of the nucleon, 
will soon be reported by experiments [501 HI] • This quan- 
tity is calculated on the lattice in much the same way as 
the DIS structure function moments are calculated. 

In this paper we report our lattice numerical calcula- 
tions of the following four moments of the structure func- 
tions: the quark momentum fraction, (x)g(//), the helic- 
ity fraction, {x)Aq{^J'), the tensor charge, and 
twist-3 coefficient df of the 32 polarized structure func- 
tion. These are the moments that can be calculated with- 
out finite momentum transfer. For the former three mo- 
ments, the momentum and helicity fractions and tensor 
charge, we restrict ourselves to the isovector flavor com- 
bination, g = M — d, as this simplifies non-perturbative 
renormalization. All three are non-perturbatively renor- 
malized and readily comparable with the corresponding 
experiments. 

The numerical calculations of these moments use the 
lattice gauge ensembles generated by the RBC and 



I 

UKQCD Collaborations with "2+1" flavors of dynami- 
cal domain- wall fermions (DWF). The good chiral and 
flavor symmetries of DWF make our calculations and 
analyses straightforward: in contrast to more conven- 
tional fermion formalisms such as staggered or Wilson, 
there is no question in defining nucleon quantum num- 
bers, nor complications arising from explicit breaking of 
chiral symmetry. This advantage is especially important 
in non-perturbatively renormalizing the results so they 
can be compared with experiment and phenomenology. 
In this paper we report results from the ensembles with 
lattice cutoff a^^ = 1.73(3) GeV IJ. We consider two 
spatial volumes, with linear size about 2.7 and 1.8 fm 
each. The strange quark mass is fixed at a value about 
12% heavier than its physical value, and the degenerate 
up and down quark mass is varied for four values from 
about three quarter to one fifth of the strange quark 
mass. Since we only vary the light quark mass in our 
simulation while the strange quark mass is held fixed, in 
the following we call the light up and down quark mass 
TO/, in lattice units, unless explicitly stated otherwise. 

The rest of this paper is organized as follows: We ex- 
plain our computational method in Sec. |Tlj In Sec. |III| 
we first summarize the numerical lattice QCD ensembles 
used for this work. Then we discuss in detail the known 
systematic errors in the relevant form factors calculated 
on these ensembles. The numerical results are presented 
in section IV Finally, we give the conclusions in Sec. |Vj 
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We note that some preliminary results from this study 
were presented in Refs. p2H24] . 



II. FORMULATION 



TABLE L Operators used in the structure function moment 
calculations, including the notation for the operator, the ex- 
plicit operator form, the hypercubic group representation, the 
correlator ratios and the projection operators used in the non- 
perturbative renormahzation of the operator in Eq. (|17[). 



We refer the reader to our recent publications 
and references cited there in for details of our computa- 
tional method. Here we give a brief summary for read- 
ers' convenience. We use the standard proton operator, 
B = £abc{UaCj5db)uc to Create and annihilate proton 
states. We Gaussian-smear this operator for better over- 
lap with the ground state with both zero and finite mo- 
mentum. A Gaussian radius of 7 lattice units was cho- 
sen after a series of pilot calculations. Since the up and 
down quark mass are degenerate in these calculations, 
isospin symmetry is exact. This is of course a well-known 
good approximation. We project onto the positive-parity 
ground state, so our proton two-point correlation func- 
tion takes the form 



1 



It 



source )), (9) 



with t — isink — ^sourcc- We insert an appropriate operator 
0(q,t') at time t' , ^source ^t' < tsink, and possibly finite 
momentum transfer q, to obtain a form factor or struc- 
ture function moment three-point correlation function. 



Cspt {t^t'^Q) — / , ^ap{Bi3{tsink)0{q,t')Ba{tsourcc))i 

(10) 



Q,/3 



with appropriate projection, F 



l + 7t 



for a spin- 



1 + 7* 

unpolarized, and F = — - — ijs'jk, k A, for a polarized 

nucleon. Ratios of these two- and three-point functions 
give plateaux for < t = t' — tsouicc < t that are fitted 
to a constant to extract the bare lattice matrix elements 
of desired observables: e.g. at = 0, we use the ratio 



(O) 



bare 



'-'3pt 



it,r) 



(11) 



In this paper we limit ourselves to those low structure 
function moments that are calculable at = 0, as are 
listed in Table U 

The structure function moments are renormalized non- 
perturbatively using the Rome-Southampton regulariza- 
tion independent (RFMOM) scheme [HI [51]. The chiral 
symmetry of DWF is relied on to suppress mixing with 
lattice-artifact operators. The specific procedures for the 
operators studied in this paper have been described in 
previous RBC publications [25||26]. Here we summarize 
them for the readers' convenience. First, we take the 
Fourier transform of the Green's function for operator 
Or constructed from a point source propagator at the 



quark momentum fraction {x)q 
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hypercubic group rep. 
correlator ratio ^(a;). 
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quark helicity fraction {x)Aq 
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transversity {l)^. 
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(V;(x)Or(0)V(y)),(12) 



{S{x,0)TS{0,y)), (13) 



where Or is one of ©14, Ogl, 0^^4j or O'^^^y The needed 
Fourier-transformed point-source and point-split-source 
propagators are 

^(p;a)=^e-^P--5(a;;0) (14) 

X 

D.Sip; a)^Y. V^'"'^ [^(^; -A)C/m(-A) 



5(x;A)C/^(0)] 



(15) 
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Next the Green's function is amputated, and evaluated 
for the case of exceptional momenta (p = p'), 

Aorb;a) = {S{p;ay^)Gor{p;p;a,){S{p;ar'^)- (16) 

Z^^ is obtained by requiring the renormalized Green's 
function, after some suitable projection, be equal to its 
tree-level counterpart pS], 

Zor{fJ.;a)~^Zq{fi;a) = ^Tr (Aor (p; a)^'r)|p2=^2 .(17) 

The projectors for each Or are listed in Table |l] Since we 
wish to match these renormalized operators to the per- 
turbative MS scheme, the renormalization scale must 
be large enough for perturbation theory to be valid, but 
not so large to introduce lattice artifacts. Thus /i should 
satisfy Aqcd ^ <C 1/a. In practice, we have found 
that the upper bound is not so strict, and can be re- 
placed by the milder condition that (pa)^ < 3. 

Finally, the following steps, similar to those from [55], 
allow us to convert the renormalization constants to the 
continuum MS scheme at 2 GeV. 

1. Obtain Z^^(/i): The ratio of Zor(p'-,o.)/Zq{iJ,;a) to 
ZA/Zq{fi; a) is computed and yields Zop(/^; a)/Z^. 
Each of the factors in the ratio is first extrapolated 
to the chiral limit, mf — — TOi-os, at fixed momen- 
tum. Using Z A ^0-7161 [J, we can determine 



Zor iP'] a) in Eq. ( 17 ) , the renormalization constant 
in the RI scheme, which we denote as Z^^. 

2. Convert to MS scheme: We are interested in con- 
tinuum quantities, mostly calculated in the MS 
scheme. The conversion factors between RI and 
MS schemes for the operators discussed here have 
been calculated in Refs. [30l |3l]. To get Z^{ii), 

we use ah (ij.) obtained by numerically solving 
the renormalization group equation with the four- 
loop anomalous dimension |32j and initial condition 

ai^\mz) = 0.1176 [33], following the method in 
Appendix A of Ref. [34] . 



3. Running to 2 GeV with two loop anomalous dimen- 
sions [31] |35| . This will take away the continuum 
running factor. 

4. Remove (ap)^ lattice artifacts: We fit the remaining 
momentum dependence, which we will observe to 
be very small, to the form / = A{ap)'^ + B and 
finally get Z^(2GeV). 



III. ENSEMBLES 

We follow the same sampling procedure as in our nu- 
cleon form factor calculations reported in Ref. [27] . 



TABLE II. A'^conf, A'^Bcp and Nmcas denote number of gauge 
configurations, trajectory separation between measurements, 
and the number of measurements on each configuration, re- 
spectively, on the 24^ ensembles. The table also lists the pion 
and nucleon mass for each ensemble |27| . 



ruf 


conf 


A^scp 


mcas 


m,[GeV] A/jv[GeV] 


0.005 


93^ 


10 


tl 


0.3294(13) 1.154(7) 


0.01 


356 


10 


4 


0.4164(12) 1.216(7) 


0.02 


98 


20 


4 


0.5550(12) 1.381(12) 


0.03 


106 


20 


4 


0.6681(15) 1.546(12) 



^ Total number of configurations is actually 646. We carry out 
extra measurements on a subset of these (286 configurations) to 
improve the statistics using different source positions. 

^ Two measurements with the double-source method gives 
effectively four measurements. 



A. Statistics 

The RBC-UKQCD joint (2-hI)-flavor dynamical DWF 
ensembles [1] are used for the calculations. These ensem- 
bles are generated with Iwasaki gauge action [33] at the 
coupling /3 — 2.13 which corresponds to the lattice cutoff 
of a^^ — 1.73(3) GeV, which is determined from the 
baryon mass [1 . 

The dynamical up, down and strange quarks are de- 
scribed by DWF actions with fifth-dimensional extent of 
Ls — 16 and the domain-wall height of M5 = 1.8. The 
strange quark mass is set at 0.04 in lattice units and 
turned out to be about 12% heavier than the physical 
strange quark, after taking into account the additive cor- 
rection of the residual mass, TOi-os — 0.00315. The degen- 
erate light quark masses in lattice units, 0.005, 0.01, 0.02 
and 0.03, correspond to pion masses of about 0.33, 0.42, 
0.56 and 0.67 GeV and nucleon masses, 1.15, 1.22, 1.38 
and 1.55 GeV. 

Two lattice sizes are used for our study, 16'^ x 32 and 
24^ X 64, corresponding to linear spatial extent of approx- 
imately 1.8 and 2.7 fm, respectively. The smaller volume 
ensembles, calculated only with the heavier three light 
quark masses, are used for a finite-volume study. On the 
16^ ensembles we use 3500 trajectories separated by 5 tra- 
jectories at m/ = O.OI and 0.02, and by 10 at 0.03. The 
main results are obtained from the larger volume ensem- 
bles with the number of the configurations summarized 
in Table HH 

On the larger volume, at the heavier three quark 
masses we make four measurements on each configura- 
tion with the conventional single source method using 
tsrc = 0, 16, 32, 48, or 8, 19, 40, 51. At the lightest mass 
the double-source method [17] is used, and two measure- 
ments on each configuration are carried out using the 
source pairs of (0, 32) and (16, 48), or (8, 40) and (19, 51). 
We made an additional two measurements on roughly 
half of the configurations with one or the other source 
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pair. This means that we make four, double-source mea- 
surements on almost half of the configurations, while two 
double-source measurements are carried out on the re- 
maining configurations. 

On the smaller volume, a single source is used, however 
the location of this source is shifted for each successive 
measurement in the order (x,y,z,t)=:(0, 0,0,0), (4,4,4,8), 
(8,8,8,16), (12,12,12,24), reducing autocorrelations. 

In order to reduce possible auto-correlations among 
measurements, they are averaged on each configuration 
and then blocked into bins of 40 trajectories for the 24"^ 
ensembles, and 20 trajectories for the 16'^ ensembles. The 
statistical errors are estimated by the jackknife method 
on the blocked measurements. 

Finally, the non-perturbative renormalization con- 
stants were computed on the four 24^^ ensembles, on 
roughly 50 configurations each, separated by 40 trajecto- 
ries. The maximum momentum value in units of in/Li 
in each direction was 6 (spatial) and 17 (temporal), such 
that {paf < 3. 
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B. Correlation functions 

The quark propagator is calculated with an anti- 
periodic boundary condition in the temporal direction, 
and periodic boundary conditions for the spatial di- 
rections. We employ gauge-invariant Gaussian smear- 
ing [37l [38] at the source with smearing parameters 
{N, uj) = (100, 7) which were chosen after a series of pilot 
calculations. For the calculation of the three-point func- 
tions, we use a time separation of 12 time slices between 
the source and sink operators to reduce effects from ex- 
cited state contributions as much as possible. 



C. Systematic errors 

There are several important sources of systematic error 
that need be considered: finite spatial size of the lattice, 
excited state contamination, and non-zero lattice spac- 
ing. A chiral-perturbation-theory-inspired analysis of the 
former for meson observables suggests the dimensionless 
product, mjri, of the calculated pion mass m,r and lat- 
tice linear spatial extent L, should be set greater than 4 
to ensure that the finite-volume correction is negligible 
(below one percent), and the available lattice calculations 
seem to support this. While our present parameters sat- 
isfy this condition, we discovered that even our larger 
volume of (2.7 fm)^ is insufficient for calculating such im- 
portant nucleon properties as the axial charge |39j and 
form factors [27]: As we reduce the light quark mass, 
eventually to = 330 MeV and mj^L ^ 4, finite-size 
effects become severe, exceeding 10 %, at least for these 
quantities that are measured in elastic processes. Similar 
finite-size effects may influence the moments of structure 
hmctions we are discussing in this paper (studies of finite 
volume effects in chiral perturbation theory can be found 



FIG. 1. A nucleon observable, isovector quark momentum 
fraction, {x)u-d, from RBC 2-flavor dynamical DWF ensem- 
ble with VTiud ~ 0.02 26 , with source-sink separation of 10 
and 12: a clear systematic difference is seen. The shorter 
source-sink separation is not manifestly free of excited-state 
contamination. 



in [40l E]). On the other hand since these moments 
are extracted from experimental observables measured 
in very different inelastic processes, the finite-size effect 
may enter differently. It is an important goal of this work 
to investigate such finite-size effects on the moments of 
structure functions. 

In order to reduce contamination from excited states, 
it is important to adjust the time separation between the 
nucleon source and sink appropriately so the resultant 
nucleon observables are free of contamination from ex- 
cited states. The separation has to be made longer as 
the quark mass is decreased. In our previous study with 
two dynamical flavors of DWF quarks [26] with a similar 
lattice cutoff of about 1.7 GeV, we saw systematic dif- 
ferences between observables calculated with the shorter 
time separation of 10, or about 1.16 fm, and longer 12, 
or 1.39 fm: the differences amount to about 20 %, or 
two standard deviations (see Fig. [T]) This would suggest 
that at the shorter time separation of about 1.2 fm, the 
excited-state contamination has not decayed sufSciently 
to guarantee correct calculations for the ground-state ob- 
servables [53]. While it is desirable to use a longer sep- 
aration, it cannot be made too long in practice without 
losing control of statistical errors. In Fig. [2] we present 
the nucleon effective mass at the lightest quark mass, 
nif — 0.005. The nucleon signal begins to decay at 
t = 12, or about 1.4 fm: this is about longest distance we 
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FIG. 2. Effective mass of the nucleon correlation function 
with Gaussian smearing applied at both source and sink, for 
quark mass ruf = 0.005. 



can choose without losing the signal. As will be shown in 
detail in this paper, the bare three-point function signals 
for this source-sink separation of < = 12 are acceptable. 
We note that recently the LHP Collaboration has also 
looked at this issue in some detail [12] and ends up using 
a shorter separation of about 1.2 fm. 

For low energy quantities like the pseudoscalar de- 
cay constants, the kaon B-parameter, and the baryon 
mass, the effect of non-zero lattice spacing was estimated 
to be less than 4% for the configuration ensemble used 
in this work fl], and subsequently confirmed on a later 
ensemble with smaller lattice spacing [43, 44|. We ex- 
pect that similar errors hold for the quantities discussed 
in this paper. 



IV. NUMERICAL RESULTS 
A. Quark momentum and helicity fractions 

Let us first discuss the ratio, {x)u-d/ i^) Au-Ad, of the 
isovector quark momentum fraction to the helicity frac- 
tion. The momentum fraction, {x)u-d, which is the first 
moment of the Fi_2 unpolarized structure functions, and 
the helicity fraction, (a;)Au-Ad: which is the first moment 
of the gi polarized structure function, share a common 
renormalization because they are related by a chiral ro- 
tation and the DWF action preserves chiral symmetry to 
a high degree. Thus, this ratio calculated on the lattice is 
naturally renornialized, much like the form factor ratio 
[27], gA/gv, and is directly comparable with the value 
obtained from experiment. 

The results of our calculation are shown in Fig. [3] 
They do not show any discernible dependence on the 
up/down quark mass, outside of the statistical error bars, 



FIG. 3. Ratio of the bare, isovector, momentum and helicity 
fractions, {x)u-d/{x}Au-Ad, which is naturally renormalized 
for DWF. Both volumes are shown, (2.7 fm)^ (circles) and 
(1.8 fm)3 (squares). The square symbols have been moved 
slightly in the plus x-direction. They are in good agreement 
with experiment which is denoted by the star. No discernible 
dependence on volume nor pion mass can be detected. 



and are in good agreement with experiment. This is 
in contrast to the renormalized ratio of gA/dv of elas- 
tic form factors which at the lightest point deviates sig- 
nificantly from heavier mass results and the experiment 
as a result of a large finite-size effect [25 . This sug- 
gests the moments of inelastic structure functions such 
as the momentum fraction, {x)u-d, and helicity fraction, 
(x) Au-Ad, may not suffer so severely from the finite-size 
effect that plagues elastic form factor calculations. In- 
deed the results obtained from the smaller (1.8 fm)'^ vol- 
ume, also shown in Fig. [3| do not deviate significantly 
from the constant behavior of the larger volume results, 
albeit with larger statistical errors. 

Next we discuss the absolute values of the isovector 
quark momentum fraction, {x)u-d- This is the first mo- 
ment of the unpolarized structure functions, Fi and F2. 
In Fig. |4| we show the bare lattice matrix elements as 
ratios of three- and two-point functions for the two light- 
est quark mass values of m/ = 0.005 (circles) and 0.01 
(squares) . We extract bare values of the desired matrix 
element by averaging over time slices 4 to 8 (values are 
summarized in Tables IV andrvl). 



These bare values need be renormalized in order to 
be compared with experiment. In Fig. [5] we present 
the non-perturbatively determined renormalization for 
the operator Of^. The filled circles are the renormaliza- 
tion constants in the RI-MOM scheme at scale /i^ — P^, 
which is not scale independent. The filled squares corre- 
spond to the renormalisation constant given in the MS 
scheme at /i = 2 GeV, where there remains only residual 
scale dependence proportional to (ap)^ lattice artifacts. 
After removing the remaining (ap)^ dependence as de- 
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FIG. 4. Signals for the ratio of three- and two-point functions 
for the bare quark momentum fraction, {x)u-d- Quark mass 
0.005 (circles) and 0.01 (squares). 
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FIG. 5. Non-perturbative renormalization factor for the 
quark momentum fraction, {x)u-d ■ Circles denote the RI- 
MOM values, squares the MS ones. The line denotes a linear 
fit used to remove the leading 0({ap)^) lattice artifacts. 



scribed in Sec. [TTj we obtain a renormalization factor of 
Zgf^(2 GeV) = 1.15(4). 

Using this renormalization constant, the quark-mass 
dependence of the momentum fraction is shown in Fig. 
[6j The results from the smaller (1.8-fm)^ volume for 
the heavier three quark mass values are in agreement 
with respective mass-value results from the larger vol- 
ume. These heavier points stay roughly the constant 
which is about 70 % higher at ~ 0.26 than the exper- 
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FIG. 6. Renormalized quark momentum fraction, {x)u-d- 
Both volumes are shown, (2.7 fm)'' (circles) and (1.8 fm)^ 
(squares). The square symbols have been moved slightly in 
the plus x-direction for clarity. 



iment, about 0.15. This behavior is not so different from 
old RBC quenched results [IS] and other recent ones [7] 
with similar up/down quark mass. 

On the other hand, the lightest point on the larger 
(2.7 fm)'^ volume, shows a sign of deviation away from 
this constant behavior. In contrast to the form factor 
deviations that move away from experiment f27J , this one 
trends toward the experimental value. Since a lighter 
quark can more easily share its momentum with other 
degrees of freedom, this trending toward the experiment 
may well be a real physical effect: It is not necessarily a 
result of the finite spatial size of the lattice. 

Indeed, it is noteworthy that the mT^L value of 3.8 
for nif = 0.01 at L = 1.8 fm is smaller than that of 
4.5 for nif ~ 0.005 at L — 2.7 fm. In other words, 
if there would be such a finite-size effect for this quark 
momentum fraction that scales with rrij^L as seen in the 
form factors, the result from m/ = 0.01 at L = 1.8 fm 
should move away from that of to/ = 0.005 at L — 2.7 fm. 
We note that in [3D] it was predicted that finite volume 
effects in {x)u-d would only become noticeable for very 
light quark masses. 

The isovector quark hclicity fraction, (x) au-AcIi ap- 
pears as a leading twist moment of the polarized struc- 
ture functions gi and 52 ■ Figure [7] presents typical bare 
signals of this quantity on the larger (2.7-fm)'^ volume, 
for the light-quark mass points nif = 0.005 and 0.01. 
Average values are extracted as for the momentum frac- 
tion, described above. Results are summarized in Tables 
|IV| for the larger volume and Table |V| for the smaller 
volume. 

The non-perturbative renormalization of (^)a^— Ad 
(for the operator O^g^j), is presented in Fig. 8 We ob- 
tain a renormalization factor of Z^^^ (2GeV) = 1.15(3) 
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FIG. 7. Signals for the ratio of three- and two-point functions 
for the bare quark helicity fraction, {x)Au-Ad- Quark mass 
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squares the MS ones. The line denotes a linear fit used to 
remove the leading 0((ap)^) lattice artifacts. 



through the same procedure described previously. This 
value agrees very well with the corresponding value for 
the momentum fraction, as guaranteed by the chiral sym- 
metry of DWF, justifying our use of bare quantities in the 
ratio, {x)u-d/ (x) Au-Ad, earlier in this section. 

With this renormalization, (x) Au-Ad can be compared 
with the experiment (see Fig. [9j) No finite volume effect 
is apparent in the data, similar to the quark momentum 
fraction. The three heavier-mass results from the smaller 
volume again agree with the respective larger-mass re- 



0.4 



0.35 - 



0.3- 



0.25- 



0.2 - 



—I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 



<X> 



Au-Ad 



experiment 
* 



• 2.7 fm 
■ 1.8 fm 



0.15 ^ ' ' ' ■ 'I 







0.1 



I 



0.2 0.3 

m'[GeV'] 



0.4 



0.5 



FIG. 9. The renormalized quark helicity fraction, (x) Au-Ad- 
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(squares). The square symbols have been moved slightly in 
the plus x-direction for clarity. 



suits, suggesting the huge finite-size effect, seen in the 
elastic form factors, that appears to scale with rriT^L, is 
not present in this moment of this deep-inelastic struc- 
ture function, at least at the quark masses considered 
here. 

Moreover, the observable exhibits very similar quark- 
mass dependence to the momentum fraction, as can be 
expected from the near constant behavior of their ratio: 
the three heavier points stay roughly the constant and 
about 70 % higher than the experimental value, and the 
lightest point shows a sign of deviation away from this 
constant behavior. This trend toward the experimental 
value may be a real physical effect. 

Here, we note that while our results of {x)^^^ are 
in agreement with rif = 2 Wilson results [45 , they 
differ significantly from the LHP mixed-action calcula- 
tions |46j . Their values are significantly lower, by about 
20 %. The main source of this discrepancy is likely 
due to the use of perturbative renormalization by the 
LHP Collaboration. In the LHP mixed-action calcula- 
tions [46], the renormalization constant is evaluated by 
Zo = {ZoIZaY"'^- X ^^™-P«rt- foj. the operator O. 
In the same manner, we use the value of [ZoIZa)'^^'^'^- 
for the operator (^{34} from Ref. [47 and ^^Q"-?™*- — 
0.7161(1) from Ref [T], and evolve it to a renormaliza- 
tion scale of 2 GeV using the two-loop anomalous di- 
mension [261 I31| . obtaining the renormalization factor 
Z(^)^^(2GeV) = 0.873(34) in the MS scheme. 

If we use this renormalization factor instead of the non- 
perturbative one described earlier, as shown in Fig. [TOj 
our results will be consistent with the LHPC results. 
The difference between the non-perturbative and per- 
turbative renormalization factors suggests a systematic 
error of about 25% should be assigned to the latter. Fur- 
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thermore, those lightest points in both quark momentum 
and helicity fractions are quite close to the experiments, 
while the non-perturbatively renormalized ones are sig- 
nificantly away from the experiments. This indicates 
that the perturbative calculation of the renormalization 
constants significantly underestimates the renormalized 
value of these particular quantities, and then exhibits an 
accidental consistency with the experiments. 

As mentioned before, it is observed that there is a 
noticeable nonlinearity in the data of both {x)u-d and 
{x) Au-Ad- These trends toward the experimental val- 
ues are easily seen in Fig. 11 where the (2-|-l)-fiavor 
and previous RBC quenched^Jl, and 2-flavor results 
are plotted together with the leading nonlinear behav- 
ior predicted in heavy baryon chiral perturbation theory 
(HBChPT) 



C 



1 - 



3£ 



1 



{x)Au-Ad = C 



n2 



In 



In 



(18) 



(19) 



Although our lightest point may be beyond the applica- 
bility of HBChPTj^ the downward trends are expected to 
develop at least in the vicinity of the physical pion mass 
point. As a simple prediction, the curve is shown using 
experimental values for the nucleon axial charge and pion 
decay constant, qa = 1.269 and Ft^ = 92.8 MeV, a chiral 
scale fi = niN = 940 MeV, and by setting the unknown 
low energy constants to zero, e(/i) — e{ii) = 0. The 
values in the chiral limit are obtained by requiring the 
curves to agree with experiment at the physical point. It 
is interesting to note that the quenched results show no 
hint of this behavior, while the 2-flavor ones are incon- 
clusive. We remind the reader that for the 2-flavor data, 
a smaller time separation between sources, iscp ~ 1-16 
fm, was used. When the separation was increased for the 
lightest mass, the momentum and helicity fractions drop, 
but with a signiflcant increase in the statistical error. 

The physical point values {x)u-d — 0.218(19) and 
{x}Au-Ad — 0.256(23), determined by simple hnear chi- 
ral extrapolation of the three lightest points, overshoot 
the experimental values by more than 2-3 standard devi- 
ations as shown in Fig. [121 other hand, nonlinear 
fit forms motivated by HBChPT, 

{x)q,Aq{nil) ^co + ci-ml + diml log(m^), (20) 

can easily accommodate all four data points and produce 
extrapolations to the physical point, in agreement with 



^ It can be observed in some particular cases like the nucleon 
axial charge qa and the nucleon root-mean-squared (rms) charge 
radius [2611271. 



the respective experimental values, albeit with large sta- 
tistical errors. Although this indicates a favorable trend, 
definitive results require simulations with several lighter 
quark masses than our lightest one. All fits are summa- 
rized in Table rvnl 



B. Transversity (Tensor charge) 

Results for the bare isovector tensor charge, {l)su-Sdj 
are presented in Fig. |13[ and in Fig. [14] we present its 
non-perturbative renormalization constant. We obtain a 
renormalization factor of Z'^(2GeV) = 0.783(3). 

Combining them we obtain the renormalized tensor 
charge as presented in Fig. [15] and summarized in Ta- 
ble |VI[ These provide a rough physical prediction which 
is still worthwhile since the experiments are yet to report 
a value. If we fit the heavy three points with a constant 
we obtain a value of about 1.10(7). Alternatively if we 
linearly extrapolate the two lightest points we would ob- 
tain about 0.7. 



C. Twist-3 moment 

Figure[T6]presents the bare lattice signals for the twist- 
3 moment, di, of the polarized structure function 172 ■ 



They are summarized in Table VIII 



We have not yet computed the renormalization con- 
stant for this quantity. The quark-mass dependence of 
the bare values are presented in Fig. [17] Our interest 
here is in whether the perturbatively obtained Wandzura- 
Wilczek relation |5T] holds. From the smallness of the 
values obtained, we conclude it does. We note that 
our results indicate that the lightest mass points deviate 
slightly from the linear trends set by the heavier points. 

In Fig. 18 values on different volumes are compared. 
di appears to be insensitive to finite volume effects, at 
least in this range of light quark masses. 



V. CONCLUSIONS 

We have presented calculations of some of the lowest 
moments of nucleon structure functions in (2-|-l)-flavor 
QCD, using domain wall fermions and the Iwasaki gauge 
action. The calculations were carried out on two volumes 
at a single lattice spacing (a^^ = 1.73 GeV) with quark 
masses that yield pion masses in the range 0.33 to 0.67 
GeV. The results are encouraging. 

The ratio of the bare quark momentum and helicity 
fractions, which is automatically renormalized, is found 
to be independent of the light quark mass through the 
range of our calculations, ^Wstrangc < rriud < fmstrange, 
and agrees with the value obtained from experiment 
within statistical error. This is in contrast to a similarly 
automatically renormalized ratio, Qa/ 9Vj of the isovector 
axialvector and vector charges that is severely distorted 
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by the finite size of the lattice [271 EH] at the lightest 
quark mass. 

This suggests the corresponding downward trend to- 
ward the experimental values, as the quark mass de- 
creases, of both the momentum, {x)u-d, and helicity, 
(x) Au-Ad, fractions is a real physical effect. Compar- 
ison of these results on two different volumes supports 
this observation. In fact, all of the moments studied here 
agree well, within statistical errors, on two different vol- 
umes, (2.7 fm)3 and (1.8 fm)^. 

In addition to the momentum and helicity frac- 
tions, the non-perturbatively renormalized tensor charge, 
{l)su-Sd, has been computed. The chiral extrapolation. 



in particular, needs to be understood before an accurate 
prediction can be made. Since upcoming experiments 
have yet to report a value, we give a rough estimate, 
from two different chiral extrapolations, that its value 
lies in the range 0.7-1.1. 

The twist-3 moment of the (?2 structure function, di , is 
also obtained. Though yet to be renormalized, its small- 
ness suggests the Wandzura-Wilczek relation holds. 

The possibility that the long sought curvature of the 
moments in the chiral regime is becoming visible in our 
results has encouraged us to start calculations at even 
smaller quark masses (m^ w 250 and 180 MeV), on an 
even larger lattice (L sa 4.5 fm). This ensemble, which 



11 



0.35 



0.3- 



0.25- 



—I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 



<X> 



u-d 



0.2- 



0.15 



^ Linear fit (3 lightest points) 
A Linear+log fit (all 4 points) 



experiment 



0.1 



0.1 







0.4 



0.35 
0.3h 
0.25- 



■ <x> 



0.15- 



0.2 0.3 
m'[GeV'] 



0.4 



0.5 



0.1 



T 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 



0-2: ^Jlf 



Au-Ad T 



expenment 



^ Linear fit (3 lightest points) 
▲ Linear+log fit (all 4 points) 



0.1 







0.2 0.3 

m'[GeV'] 



0.4 



0.5 



FIG. 12. Linear and leading HBChPT fits to the non-perturbatively renormalized quark momentum fraction, and 
helicity fraction, (x) Au-Ad- The three lightest data points are used in the former, and all four in the latter. The diamond and 
triangle denote extrapolated values at the physical point. 



2| r 



LS- 



I- 



IS- 



—I 1 r- 



* i i i i * 



m=0.005 
m=0.01 



12 



t-t 



FIG. 13. Signals for the ratio of three- and two-point functions 
for the bare quark transversity, {l)su-Sd- Quark mass 0.005 
(circles) and 0.01 (squares). 



is being generated by the RBC and UKQCD collabora- 
tions 52J, was conceived, in part, to attain these goals 
for nucleon matrix elements. 
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TABLE III. Summary of renormalization factors in the MS 
scheme at 2 GeV in the chiral limit. 



•-^34 



-m^cs 1.149(45) 1.154(35) 0.783(3) 



TABLE IV. Bare quark momentum and helicity fractions and 
their naturally renormalized ratio on the (2.7 fm)^ ensemble. 



0.005 0.201(9) 0.240(13) 0.835(46) 

0.01 0.219(9) 0.261(14) 0.842(42) 

0.02 0.234(8) 0.286(11) 0.821(40) 

0.03 0.231(7) 0.285(10) 0.807(32) 



TABLE V. Bare quark momentum and helicity fractions and 
their naturally renormalized ratio on the (1.8 fm)^ ensemble. 

rrlf {x)u-d {x}Au-Ad {x)u-d/{x)Au-Ad 

0.01 0.221(18) 0.263(29) 0.808(89) 
0.02 0.256(14) 0.291(22) 0.875(56) 
0.03 0.236(7) 0.300(11) 0.784(27) 



TABLE VI. Isovector combination {u — d) of the quark mo- 
mentum fraction {x)q, helicity fraction {x)Aq and transversity 
{l)sq, non-perturbatively renormalized in the MS scheme at 
2 GeV. 



rrif 


{x)^s^(2GeV) 


(^)L'-Ad(2GeV) 




0.005 


0.231(14) 


0.277(17) 


1.265(45) 


0.01 


0.252(14) 


0.302(19) 


1.438(41) 


0.02 


0.269(14) 


0.330(16) 


1.384(38) 


0.03 


0.266(13) 


0.329(15) 


1.446(32) 



TABLE VII. Summary of extrapolations to the physical point 
of the renormalized first moment of the quark momentum and 
helicity fractions. 







{x)^l,{2GeY) 


{x)^LAdi2GeY) 


T • 2 

Linear vs m„ 


(3 points) 


0.218(19) 


0.256(23) 


ChPT vs ml 


(4 points) 


0.175(51) 


0.205(59) 


Experiment 




0.154(3) 


0.196(4) 
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